Testing Gravity Against Early Time Integrated Sachs- Wolfe Effect 
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A generic prediction of general relativity is that the cosmological linear density growth factor 
D is scale independent. But in general, modified gravities do not preserve this signature. A scale 
dependent D can cause time variation in gravitational potential at high redshifts and provides a 
new cosmological test of gravity, through early time integrated Sachs-Wolfe (ISW) effect-large scale 
structure (LSS) cross correlation. We demonstrate the power of this test for a class of f(R) gravity, 
with the form f(R) = — XiHq exp(— i?/A2-ffo )■ Such f(R) gravity, even with degenerate expansion 
history to ACDM, can produce detectable ISW effect at z > 3 and I > 20. Null-detection of such 
effect would constrain A2 to be A2 > 1000 at > 95% confidence level. On the other hand, robust 
detection of ISW-LSS cross correlation at high 2 will severely challenge general relativity. 

PACS numbers: 98.65.Dx,95.30.Sf 
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Introduction. — Cosmological observations provide 
unique tools to study gravity at > Mpc scales. Gen- 
eral relativity, with the aid of the cosmological constant, 
or dark energy with equation of state w ~ — 1, success- 
fully reproduces the accelerated expansion of the Uni- 
verse, indicated by SN la observations Q, along with 
the flatness of the Universe measured by the cosmic mi- 
crowave background (CMB)@| and distance measured by 
the baryon oscillations Q. However, these observational 
evidences mainly constrain the mean expansion history 
of the Universe and can be reproduced by modified grav- 
ity such as brane world DGP theory £| and generalized 
f(R) gravity^. Essentially, the large scale structure 
(LSS) of the universe, such as weak gravitational lensing, 
galaxy clustering and the integrated Sachs- Wolfe (ISW) 
effect [a, 0, M, HI liil liiJ is required to break this degen- 
eracy. 

General relativity imprints a unique signature in the 
LSS, which is scale independent linear density growth fac- 
tor D at sub-horizon scale after matter-radiation equality 
epochfi^. Modifications to general relativity not only 
changes the amplitude of D, but in general, causes D to 
be scale dependent. This unique feature of modified grav- 
ity has already been noticed in phenomenological theory 
of modified Newtonian potential[6|, |9l [lfj. It can be de- 
tected by weak gravitational lensing|f|, galaxy clustering 
0, 0| and late time ISW effect. Counter-intuitively, in 
this paper, we show that modified gravity can produce a 
detectable early time ISW effect. 

We investigate a class of f(R) gravity with action 



and field equation 
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where /r = df/dR. We design f(R) = 
— X\Hq exp(— R/\ 2 Hq), where Ai,2 are two positive di- 
mensionless constants and Hq is the Hubble constant at 
present. To mimic a ACDM universe, Ai ~ 1 is required. 
To reduce to the general relativity in the solar system 
and pass the solar system test, / -C R is required. In 
this limit, we have R — > 87rGp so i ar , where p S oiar is the 
local density where solar system tests are carried out. In 
this limit, f(R)/R~ [p c / p soial ] exp[-3p so i ar /A 2 p c ]. Since 
3psoiar/Pc ^ 10 6 ll^KPc is the critical density of the Uni- 
verse), models with A2 >C 10 6 can survive all solar system 
tests. For example, For A2 = 10 3 , this correction is of the 
order ~ 1CP 400 . Given such tiny f(R), we expect that /, 
fit, a fn and Jr-^-u in Eq. |21can all be safely neglected 
for any physical purpose. 

For the f{R) gravity, the application of Birkhoff the- 
orem to perturbations of a spherically symmetric region 
leads to scale independent D 0. We reinvestigate this 
issue by solving the structure evolution of the fully co- 
variant f(R) gravity to linear order in the metric per- 
turbation. We find that D shows nontrivial scale depen- 
dence, consistent with the results based on the Palatini 
approach [l7| . 

The H-z relation of the f(R) gravity. — Cosmological 
observations prohibit strong deviation of f(R) from the 
cosmological constant. At the limit that R(z = 0) <C 
X2H 2 , the H-z relation of f(R) gravity can have the 
same asymptotic behavior as that of ACDM. At low red- 
shift where R(a) <C X 2 H^, f(R) behaves as a cosmo- 
logical constant and the H-z relation resembles that of 
ACDM. At high redshifts where R > X 2 H§, 
and H(z) -> nj /2 (l + zfl 2 . Deviation from ACDM hap- 



pens at some intermediate redshifts where R(a) ~ X 2 Hq 
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The large scale structure of the f(R) gravity. — We will 
show that, even with this degeneracy in H-z relation and 
solar system behavior, the LSS of the f(R) gravity could 
be significantly different to that of ACDM. We choose 
the Newtonian gauge 
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FIG. 1: The H(z)-z relation and structure growth in the ex- 
ponential f(R) gravity. Top left panel: H-z. X2 —* 00 corre- 
sponds to ACDM cosmology. Top right panel: Q(k,a) cc k 2 , 
which describes the main effect of f(R) gravity to structure 
formation. We plot the result of k = 0.01/i/Mpc. Bottom left 
panel: fit (a), which determines the effective Newton's con- 
stant Gcft — G/(l + fn). For A2 ^ 100, its effect to structure 
formation can be neglected. Bottom right panel: D(k,a)/a 
(A2 = 1000), where the linear density growth factor D is nor- 
malized such that D — » a when a — > 0. 



and vanishes toward both higher and lower z. We quan- 
tify their difference by solving Eq. [21 of a flat universe to 
zero order 
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Here, a = 1/(1 + z) is the scale factor. This equation 
can be rewritten as y = Sl — C(y(a)), where y = a 3 H 2 , 
C(y(aj) = [f/6 — afn/a + Hf R ]a 3 and Q is the dimen- 
sionless matter density at present. Since C(y(a)) is com- 
pletely determined once y as a function of a is given, 
Eq. |2| can be solved iteratively by the iteration relation 
— fl - To mimic a ACDM universe, we 

fix Ai by requiring f(R(a = 1)) = -6H§(1 - Q ). The 
iteration converges quickly by taking the initial guess 
= Q + (1 - n )a 3 . For A 2 > 100, is accu- 
rate to ~ 1%. As expected, for A2 > 100, the H(z)-z 
relation is almost identical to the corresponding ACDM 
cosmology (Fig. Such f(R) gravity can not be dis- 
tinguished from ACDM by inflation, big bang nucleosyn- 
thesis (BBN), primary CMB, SN las and other measures 
of H-z relation. 



ds 2 = -(1 + 2^)dt 2 + a 2 (l + 2(t>) Y^( dxi f 
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There are four perturbation variables (j>, ip, the matter 
over-density 5 and the (comoving) peculiar velocity con- 
vergence 9. 

In general relativity, <f> = —tp, as long as there is no 
anisotropic stress. But in modified gravity, this relation 
breaks in general, ij (i ^ j) component of Eq. [2 pro- 
vides the relation between <j> and ip. For f(R) gravity, 
due to non-vanishing fR-i-j (i ^ j), 4>-il> relation becomes 
scale dependent. Throughout this paper, we neglect time 
derivative terms with respect to spatial derivative terms 
of corresponding variables. This simplification holds at 
scales k > aH < 10 _3 /i/Mpc. Since we will focus on 
the ISW effect at I > 20 and z > 3 where the relevant 
fc>5x 10~ 3 /i/Mpc, this simplification is sufficiently ac- 
curate. We then obtain 
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In Fourier space, this reads ijj = — <f>{\ — 2Q)/(1 
where Q(k,a) = -2f RR c 2 k 2 /(l + f R )a 2 and f RR = 
d 2 f/dR 2 . For clarity, we explicitly show the speed of 
light c. We will see that this scale dependent ip-<fi rela- 
tion has profound effect on the LSS. Combining Eq. [S] 
and the tt component of Eq. we obtain the new Pois- 
son equation 
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The energy-momentum tensor is still conserved and pro- 
vides the remaining two equations: 



= 0, 6 + 2H6+\v 2 il) = 
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Combining all 4 equations, we obtain the main equation 
of this paper: 
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where ' = d/da. In general relativity, Q = 0, so the 
linear density growth factor D oc S(a)/S(a,i) is scale in- 
dependent at scales k > aH/c, no matter what the form 
of dark energy is. Here, ai is the scale factor at some 
early epoch and we normalze D such that D — > a when 
a — > 0. But in f(R), the scale dependent Q(k,a) induces 
nontrivial scale dependence to D. This behavior can not 
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in ip — cj) causes a fractional CMB temperature variation 
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FIG. 2: The ISW effect. A 2 = fOOO is adopted. Top left 
panel: D/a — dD/da, which determines the sign and ampli- 
tude of the ISW effect. D is normalized such that D —* a 
when a — *• 0. Bottom left panel: the ISW effect. Bottom 
right panel: Cumulative S/N of the ISW-LSS cross correla- 
tion measurements. 



be obtained by a simple change in the effective Newton's 
constant. Furthermore, the correction Q has a nontriv- 
ial dependence on a. This is hard to realize by simply 
changing the form of the Newtonian potential (e.g. to 
Yukawa potential). 

Since Jrr < 0, there exist one apparent singularity 
Q = 2/3 in Eq. 03 where only 5 — solution is accepted 
and two at Q = 1/2,1 in the ip-(j) relation, where only 
■ijj = <p = solution is accepted. We leave this issue 
alone until the discussion section. For the moment, we 
take a modest goal by only using regions where Q < 1/2 
to constrain f(R). For A2 = 1000, this constrains us to 
region where k < 0.012/i/Mpc. 

Hereafter, we fix A 2 = 1000. At z > 1, H oc a~ 3/2 , 
D oc a 1 "" when rj = 3Q/5(2 - 3Q) < 1. Thus gravita- 
tional potential decays at high redshifts with rate oc a~ ri 
and causes an observable integrated Sachs- Wolfe (ISW) 
effect. At later time when R < \ 2 H%, Q -> (Fig. [TJ, 
the evolution of D approaches that of ACDM. For the 
exponential f(R), Q(a) peaks at z ^> 1 (Fig. [IJ, so the 
resulting ISW effect peaks at z 3> 1, as contrast to that of 
ACDM cosmology or dark energy models with w ~ — 1. 
This provides us a unique way to test this form of f(R). 
We solve Eq. [S] numerically. Initial condition is set to 
normalize D —t a when a —t 0. 

The integrated Sachs-Wolfe effect. — Time variation 
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Here, x is the comoving angular diameter distance. Since 
both ip — (j> and the LSS trace the underlying matter 
distribution, there exists an ISW-LSS cross correlation, 
with power spectrum 
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Here, <5lss is the density fluctuation of the LSS trac- 
ers, Wlss is the corresponding weighting function 
and A? • •„ r is the corresponding 3D power spec- 

(i/>-0)<5lss 

trum(variance). The above formula adopts the Limber's 
approximation, which is sufficiently accurate to serve for 
our interest at I > 20. The amplitude and sign of the 
ISW effect is determined by A\sw = D/a — dD/da. Pos- 
itive Aigw means positive correlation between ISW and 
LSS. For k > 0.007/i/Mpc, Ai SW has a bump at z ~ 6, 
whose amplitude increases towards small scales (large k) . 
This boosts early time small scale ISW signal (Fig. |3J. 

The S/N of the ISW-LSS cross correlation measure- 
ment of each I mode is 
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Here, Ccmb, Cisw ; C g are the power spectra of primary 
CMB, ISW, and galaxies, respectively, while C^b an< ^ 
(jshot are p 0wer spectra of associated shot noises, re- 
spectively. Since the exponential f(R) does not affect 
physics at z > 100, we adopt the same primordial power 
spectrum wi th p ower index n = 1, the same transfer func- 
tion BBKS [19| and the same amplitude at a; = 0.01, as 
that of the ACDM cosmology. The LSS tracers we choose 
are 21cm emitting galaxies at 3 < z < 5, which will be 
measured by proposed 21cm experiments such as Square 
Kilometer Array |20j . Singularities presented in the per- 
turbation equations limit us to I < 60, where one can 
neglect shot noises of CMB. For the estimations of LSS 
clustering signal and shot noise, biggest uncertainties are 
(1) HI (neutral hydrogen) mass function at 3 < z < 5, (2) 
21cm emitting galaxy bias and (3) specifications of 21cm 
experiments. If one adopts HI mass functions calibrated 
against observations of damped Lyman-a systems and 
Lyman limit systems, SKA can detect > 10 9 galaxies at 
z > 3 in five years across the whole sky, for a field of 
view > 10 deg 2 at ~ 300 Mhz (for details of the calcu- 
lation, see, e.g. H]]]). Detection thresholds of HI mass 
at z > 3 are > 10 9 M Q , so detected galaxies are likely 
having biases bigger than one. Then, one can neglect the 
shot noise term C^ hot with respect to C g . Taking the 
fact that C CMB > C ISW (Fig. EJ, the S/N of each I is 
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simplified to 

(S\\_ (2^ + l)/ sky r 2 

\N J ~ Cp MB /Cp w ' 1 J 

Here r is the cross correlation coefficient between ISW 
and LSS. Since r has very weak dependence on galaxy 
bias, the estimation presented here is weakly model de- 
pendent. We disregard signals from I < 20, to reduce 
confusions of ACDM cosmology or dark energy models. 
For sparse galaxy sampling which is sufficient for our pur- 
pose, SKA is able to cover the whole sky. So we assume 
that / s ky = 1. The cumulative Ylw^i^/N) 2 ^ s s h° wn m 
Fig0 

The ISW signal peaks at z > 3 and increases toward 
high I. This is hard to mimic by ACDM, dark energy 
or many forms of modified gravity. (1) For ACDM or 
dark energy models with w < — 1, at z > 3, the ISW 
effect effectively vanishes. Fig. [21 shows that ACDM 
can be distinguished from the A 2 = 1000 f(R) grav- 
ity with > 2cr confidence by the ISW-21cm emitting 
galaxy cross correlation. (2) For dark energy models 
with w > — 1, Aisw does not decrease as fast as that 
of ACDM. But the ISW signal (including contributions 
from dark energy fluctuations) decreases toward high I 
[22]| and one does not expect a detectable ISW effect. 
(3) DGP preserves the property of scale independent D 
H, 1^3, so the ISW signal decreases toward high i, like 
the dark energy case. Therefore we do not expect a de- 
tectable signal at I > 20 and z > 3. (4) For generalized 
/ oc (aR 2 + PR ab R ab + lRabcdR abcd )- n (n > 0), the 
ISW effect vanishes at high z because the / correction 
decreases much faster than the exponential f{R). So we 
expect that null detection of ISW-LSS cross correlation 
at Z > 20 and z > 3 would constrain A2 to A2 > 1000 at 
> 2(7 confidence level. On the other hand, a detection of 
such cross correlation would present as a severe challenge 
to general relativity. 

Discussion. — The scale dependence of D, as an unam- 
biguous signature of modified gravity, can in principle be 
measured from weak gravitational lensing by the mean 
of lensing tomography. Since <fi is no longer equal to — tp, 
we provide the general form of the lensing transformation 
matrix Aij 

An - % = d X (4 - 1>hjW(x, Xs) , (13) 
Jo 

where W(\iXs) — xO- ~ x/Xs) is the usual lensing ker- 
nel. All basic lensing theorems remain unchanged. For 
example, lensing shear field is still curl free [if neglecting 
second order corrections such as Born correction). For 
f(R) gravity, relation between the lensing convergence 
k = 1 — (An + A 2 2)/2 and the matter over-density re- 
sembles that of the general relativity, with 

k = ^ 2 r» f 5a- x W( X , X .)(l + In^dx ■ (14) 



It is interesting to see how well weak lensing alone can 
constrain modified gravity. For the exponential f(R), one 
complexity is that lensing mainly probes LSS at z < 1, 
where Q is small and the deviation from a scale indepen- 
dent D is small, so the constraints may be weak. This can 
be significantly improved by gravitational potential re- 
constructed from primary CMB. Combining lensing and 
CMB measurements, it is very promising to measure the 
evolution of the gravitational potential between z = 1100 
and z ~ robustly. This will put strong constraints on 
the nature of gravity. Unfortunately, due to singularities 
in the perturbation equations, we are limited to scales 
k < 0.012/i/Mpc or I < 20 at z < 1 (for A 2 = 1000). In- 
formation contained in this region is very limited and 
could be contaminated by other physics such as dark 
energy fluctuations. Solving the field equation crossing 
those singularities consistently is nontrivial. We leave 
this work for future study. 

The Q = 1/2,2/3,1 singularities may be caused by 
awkward gauge choice, the neglecting of time derivative 
terms with respect to corresponding spatial derivative 
terms, or the failure of the perturbation approach. For 
example, for Q = 2/3, the only solution <5 = does not 
depend on initial conditions. This could be caused by 
neglecting time derivative terms, which erases some de- 
grees of freedom. These issues require detailed study. 
But if these singularities in LSS equations are physical, 
they can be applied to rule out many forms of modified 
gravities as alternatives to dark energy or general rela- 
tivity. To produce a similar expansion history as those 
of dark energy model, (1) R should increase when a de- 
creases and (2) f(R(a = 1)) should be negative in order 
to mimic positive dark energy. Furthermore, in order not 
to affect inflation, BBN and primary CMB, /(#(a -> 0)) 
must be sufficiently small. A sufficient (but not nec- 
essary) condition satisfying the BBN constraint is that 
f(R(a -> 0)) -> 0. The exponential f(R) and 1/R n 
f(R) all fall into this class. This results in fn > at 
least at some early epoch a+. As we have seen from 
previous discussions, fjm < is a sufficient condition 
for the existence of singularities. To avoid singularities, 
fun > must be satisfied at all epochs. However, we will 
see that this requirement contradicts with (1) and (2). 
fun > results in f R (a < a + ) > fu(a + ) > 0, because 
R(a < a_|_) > R(a + ). So, / increases toward high red- 
shift, crosses over zero at some epoch and then increases 
more quickly than R. Since when a — > 0, R cx a~ 3 , / 
increases more quickly than a -3 and thus more quickly 
than the matter density. This violates condition (2). It 
could have non-negligible effect on BBN and contradicts 
our expectation. On the other hand, only for those f(R) 
gravities in which f(R(a — > 0)) does not vanish, singu- 
larities in LSS equations can be avoided. A \ogR f(R) 
gravity is such an example. 

To demonstrate the power of LSS to constrain grav- 
ity, we adopt a conservative requirement to avoid sin- 
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gularities at k < k s . At the limit that A2 1, Q 
peaks at a = (2A2/9^o) ^ 3 and the peak amplitude is 
~ 12(1 - n )(2/9n e) 2 / 3 \~ i/3 (ck/ H ) 2 , where we show 
the speed of light c explicitly. To avoid singularities at 

A ^ 2 - 5 * ios yy 3/2 <i5) 

should be satisfied. 
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